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ABSTRACT 

Wc report the first extensive set of optical photometric observations of the counterpart 
to SAX J1808. 4-3658 (V4580 Sagittarius) in quiescence. The source was detected at 
V ^ 21, 5 magnitudes fainter than at the peak of its 1998 outburst. However, a 
comparable ~6% semi-amplitude 2hr modulation of its flux is revealed. This has the 
same phasing and approximately sinusoidal modulation as seen during outburst, and 
with photometric minimum when the pulsar is behind the companion. The lack of a 
double-humped morphology rules out an ellipsoidal origin, implying that the bulk of 
the optical flux does not arise from the companion. Moreover, applying crude modelling 
to the disc and X-ray irradiated face of the donor shows that the internal energy release 
of a remnant disc (with mass transfer driven by gravitational radiation) is sufficient 
to explain most of the optical emission, and with the modulation due to the varying 
contribution of the heated star's face. We note that this model is also consistent with 
the much lower X-ray to optical flux ratio in quiescence versus outburst, and with the 
phasing of the optical modulation. 

Key words: accretion, accretion disks - binaries: close - pulsars: individual: SAX 
J1808. 4-3658 - stars: individual: V4580 Sagittarius - stras: low-mass, brown dwarfs 
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1 INTRODUCTION 

The low-mass X-ray binary (LXMB) transient SAX 
J1808. 4-3658 was first detected (in outburst) by the 
5eppoSAX Wide field Camera in 1996 September 
(In't Zand et al. 1998), following a non-detection in 
a preceding August observation. Three type-I X-ray 
bursts were also observed demonstrating the presence 
of a neutron star accretor. The third exhibited clear 
double-peaked morphology implying an Eddington- 
limited photospheric radius expansion event, yielding 
a distance estimate of 2.5kpc (In't Zand et al. 2000). 

During another outburst in 1998 April/May, SAX 
J1808.4-3658 was detected by RXTE, enabling the mea- 
surement of the fastest ever LMXB X-ray pulsation, 
with a period of 2.49ms (Wijnands & van dcrKlis 1998). 
Chakrabarty & Morgan (1998) utilized precise tim- 
ing of the millisecond pulses to determine many pa- 
rameters of its binary orbit: period, epoch of mean 
longitude, projected semi-major axis, eccentricity and 
pulsar mass function. Hence, they were immediately 



able to place constraints on the nature of the sys- 
tem, in particular the secondary's spectral type, sug- 
gesting a very low mass (< 0.1 Mq), but irradiation- 
bloated/ablated, dwarf star. Lastly, Roche et al. 1998 
were able to identify an optical counterpart consis- 
tent with the -SeppoSAX X-ray position, even within 
the smaller error circle yielded by the RXTE X-ray 
pulse timing. We note that this counterpart has recently 
been designated as the variable star V4580 Sagittarius 
(Kazarovcts et al. 2000). 

Additional constraints on the secondary can be 
placed by observations of the optical counterpart. Giles, 
Hill & Greenhill (1999) report on CCD photometry per- 
formed with the Mt Canopus Im telescope during de- 
cline from the 1998 outburst in April-June. Their cov- 
erage consists mostly of single pointings a few nights 
apart, in order to follow the long-term lightcurve, which 
showed the source to fade from V=16.72 to around 18.5, 
1 and 3 weeks after the start of the X-ray decline. A 
later observation on June 27 failed to detect the source. 
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corresponding to V >20.5. They also obtained short 
l/-band sequences on three nights, but only amounting 
to ^0.25, ~1.0 and '^1.2 cycles of the 2 h period re- 
spectively. From these sparse data they concluded that 
there was evidence for a 0.06 mag semi-amplitude mod- 
ulation at the 2 h X-ray period, with phasing such that 
the minimum occurs when the companion lies between 
the observer and pulsar. 

In order to place further limits on the secondary, we 
obtained the first high time-resolution CCD photome- 
try of the counterpart in quiescence, using the 1.9m 
telescope at the South African Astronomical Observa- 
tory (SAAO), and it is these results that we present in 
this paper. 



2 OBSERVATIONS AND DATA 
REDUCTION 



N 




20' 



Figure 1: Finding chart for V4580 Sgr (SAX J1808.4- 
3658) in quiescence. This is a co-add of 13 B-band im- 
ages and totals 52 min of exposure. Note the faintness of 
the target and and the closeness of neighbouring stars. 
The three local standard stars used for calibration are 
indicated. 

Observations of a small (50 x 33 arcsecs) region 
surrounding the optical counterpart to SAX J1808.4- 
3658 were made using the UCT-CCD fast photome- 
ter (O'Donoghue 1995), at the Cassegrain focus of the 
1.9m telescope at SAAO, Sutherland on 1999 August 10 
and 2000 July 3. The UCT-CCD fast photometer is a 
Wright Camera 576x420 coated GEC CCD which was 
used here half-masked so as to operate in frame trans- 
fer mode, allowing exposures of as short as 2s with no 
dead-time. However, given the faintness of the quies- 
cent counterpart we employed much longer integration 
times. In 1999, we concentrated on obtaining a high 
signal-to-noise lightcurve, and therefore we opted for 
a continuous time-series of 30s white light (no filter) 



exposures. In 2000, we were again favoured with good 
seeing and stable photometric conditions initially, so we 
attempted to obtain colour information on the binary 
modulation. In fig. |l| we show a deep (52 min expo- 
sure) B band finding chart, constructed from the 13 
shifted and co-added images from the interval of best 
seeing during these observations. This demonstrates the 
difficulties in obtaining good photometry of this faint 
target, as it is quite severely crowded; only with good 
seeing, < 1" (indeed truly excellent for this telescope) is 
there any possibility of resolving the counterpart from 
its near neighbour. However, the seeing did deteriorate 
and become more variable towards the end of the run. 
In this instance, we cycled through BVR filters, with 
typical exposure times of 240, 180, 150s in each respec- 
tively. For the 2 h period, this still yields approximately 
0.1 phase resolution. 

We performed data reduction using IRAF, includ- 
ing photometry with the implementation of DAOPHOT 
II (Stetson 1987). Point spread function (PSF) fitting 
was always employed, as it is essential for obtaining 
the best possible photometry, in such crowded condi- 
tions. The details of this procedure are given in Homer, 
Charles & O'Donoghue (1998). Lastly, in order to re- 
duce systematic effects differential photometry was ap- 
plied, wherein the magnitudes of a star are calculated 
relative to an ensemble of local standards (bright, non- 
variable stars). 

3 RESULTS 

3.1 Colour photometry 

For the 1999 August 10 white light observations no 
strict calibration is possible. However, if we assume that 
the unfiltered passband of the CCD corresponds roughly 
to V , we may use observations of the E-region standard 
E752 (Menzies & Laing 1980) to estimate that the coun- 
terpart had 20 on that date, consistent with the 
limit of Giles, Hill & Greenhill (1999) given the uncer- 
tainties. 

For the sets of colour images from 2000, we use 
the standard star E508 (chosen for its similar colours, 
and comparable airmass when observed), and aperture 
photometry to derive the calibrated magnitudes of lo- 
cal standard stars within the small field of the CCD 
(see fig. |l|). Finally, the magnitude of the target is cal- 
culated from its differential magnitude with respect to 
these local standards, these measurements necessarily 
making use of PSF fitting. We find B = 22.0 ± 0.1, 
V = 21.5±0.1 and R = 20.9±0.1 all uncorrected for red- 
dening. The magnitude errors quoted arise largely from 
counting statistics as the target is simply so faint, to- 
gether with uncertainties in the transformation between 
aperture and PSF magnitudes. Compared to its colours 
during outburst (Roche et al. 1998) of 16.6 ± 0.2 
and R = 16.1 ± 0.2, i.e. {V - R)hi = 0.5 ± 0.3 cf. 
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Tabic 1: Quiescent optical properties of V4580 Sagittar- 
ius/SAX J1808.4-3658 



Quantity 




Lower limit 


Upper limit 


Apparent mag. 




17V 
11 . 1 


zl.U 


(de-reddened) 


Vo 


18.2 


20.8 




Rq 


17.7 


21.0 


Absolute 




5.7 


9.0 


magnitude 


My, 


6.2 


8.8 




Mk, 


6.3 


8.3 


Flux density 


Fbo 


2.5x10"" 


5.3x10-1*^ 


( erg cm~^ s^^ 


Fvo 


1.8x10-1^ 


2.0x10-1*^ 






1.3x10-1^ 


8.4x10-1^ 



{V — R)io ~ 0.6 ± 0.2; apart from being 5 magnitudes 
fainter, they are consistent within the uncertainties, 
agreeing with the unchanging V — I values found by 
Giles, HiU & GreenhiU (1999) during the decUne itself. 

Furthermore, in order to facilitate other compar- 
isons we may attempt to correct for interstellar ex- 
tinction towards SAX J1808. 4-3658. We utiUze the col- 
umn density estimated from RXTE X-ray spectral fits 
by Heindl & Smith (1998) to set an upper limit, and 
the galactic value in the direction of SAX J1808. 4-3658 
for the lower limit (Dickey & Lockman 1990), to de- 
rive 6.3 > Nh > 1-3 X 10^^ cm"^. Applying the rela- 
tions Nh = 1.79 X lO^i^y (Predehl & Schmitt 1995), 
E{B-V) = Ay/3.2, Ab = 1.32Ay, and An. = 0.81Ay, 
we may estimate corrected colours of —0.6 < [B—V)q < 
0.3 and -2.4 < (V - R)o < 0.5. The corrected flux 
densities can also be obtained using the calibration for 
an AOV star in Drilling & Landolt (2000), as may the 
absolute magnitudes assuming a 2.5 ± O.lkpc distance 
(In't Zand et al. 2000). All these data are summarised 
in table ^ 

Lastly, independent of both distance and redden- 
ing estimates we may compare the X-ray to optical 
flux ratios in outburst and quiescence. We use the out- 
burst RXTE/PCA X-ray and contemporaneous opti- 
cal observations on 1998 April 18 (Gilfanov et al. 1998; 
Roche et al. 1998), whilst we adopt an intermediate 
value between the two quiescent X-ray measurements by 
BeppoSAX on 1999 March 17-19 (Stella et al. 2000) and 
ASCA on 1999 September 17 (Dotani, Asai & Wijnands 
2000) together with our own optical data of 2000 July. 
Hence, the quiescent value will be uncertain by a factor 
of a few, given the known long-term X-ray flux vari- 
ations. The X-ray fluxes have been transformed to the 
2-10 keV range using W3PIMMS at HEASARC, and as- 
suming a power-law spectrum with photon-index F = 2. 
We calculate Fx (2 - 10keV)/Fy = 9.0 x IO-178.5 x 
10-13 = 1100 in outburst and 1x10-13/9 3x10-15 ^ 10 

in quiescence where Fy = 990Ax/y, the flux density 
in V. The ratio in outburst is comparable to other out- 
bursting transients and persistent LMXB, where the 



dominant optical emission is presumed to arise from 
X-ray reprocessing - clearly this is not the case in qui- 
escence. 

3.2 Orbital modulation 

-0.15 p . , . . , . . . . , . . . ^ , . . . ^ 

SAX J1808, 4-3658 (SAAO White light) 

1999 August 10 

-0,1 - 




0.1 - 



0.5 1 1.5 

orbital phase 

Figure 3: Phase folded and binned lightcurve (white 
Hght) of V4580 Sgr (SAX J1808.4-3658), with a con- 
strained sinusoidal fit over-plotted. Two orbital cycles 
are shown for clarity. 

The lightcurve from 1999 August spans a full 2 cy- 
cles, with another half cycle, where the effects of air- 
mass (and deteriorating seeing) become a slight prob- 
lem. Clearly, this is insufficient to make any independent 
measurement of the period, but we can easily check for 
consistency with the pulse timing result. In figure ^ we 
plot the linearly detrended lightcurve, along with a fit- 
ted sinusoid with Px = 2.01h and the projected X-ray 
ephemcris. Clearly, there is no doubt from our dataset 
that the quiescent counterpart exhibits a smooth ap- 
proximately sinusoidal modulation at the 2 h period, 
and the fit gives the semi-amplitude as 0.063 ± 0.003 
mag. For this fit and later analysis we opted to omit 
the last half cycle of data, where the source appears to 
reach its brightest level, as this may be due to the poorer 
seeing. We note that the amplitude, like the colours, is 
consistent with that of the high state source. The scatter 
seen in the similar lightcurves of Giles, Hill & GreenhiU 
(1999) could be due to an intrinsic property of the high 
state source (it is known on occasion to exhibit com- 
plex X-ray flaring behaviour; van derKlis et al. 2000; 
Wijnands et al. 2000), which is likely to be present in 
the reprocessed optical emission. In order to study the 
morphology more closely, we have also folded and phase 
binned the data, which is presented in figure^. The pure 
sinusoid appears to be a very good fit, only at around 
0.9 in phase does the data scatter, which corresponds 
to a region in phase for which only one cycle of data 
is present (owing to a gap in the first cycle). Quantita- 
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Figure 2: Lightcurve of V4580 Sgr (SAX J1808. 4-3658), taken in white light (fiUed circles), with a star of comparable 
brightness plotted offset below (squares). The data has been binned up a factor of two to 60s time resolution for 
clarity. Over plotted is the sinusoidal fit to the data (solid line), with only the period constrained to the orbital 
value. The dotted line indicates the X-ray ephemeris, i.e. the mean orbital longitude of the system as derived from 
the X-ray pulse arrival time delays (Chakrabarty & Morgan 1998). As during the outburst observations of Giles, 
Hill & Grecnhill (1999), the time of photometric minimum is consistent with a mean orbital longitude of 90°, that 
is when the pulsar lies behind the companion. 



tively, we obtain Xv = 1-1 (for 12 d.o.f) for the fit to 
the folded and binned data. 

In July 2000 we attempted to investigate the colour 
dependence of the two hour modulation. However, the 
target turned out to be even fainter than in August 
1999, and hence the signal-to-noise of our data is 
marginal. In fact, we found it necessary to firstly select 
the best seeing images (based on check star behaviour) , 
and then phase fold the data. The data have also been 
detrended using a quadratic polynomial constrained to 
approximate the airmass variations. In figure ^ we show 
the resulting BVR lightcurves for V4580 Sgr and a star 
of comparable brightness. Clearly, the phase folding re- 
veals apparent modulations in each band, although in 
the case of V the check star also shows only slightly 
smaller amplitude variations. Given the data quality we 
must be cautious with regard to the significance of these 
modulations. However, statistically we do find that a 
(period and phase constrained) sinusoid fit is always an 



improvement over a constant value at > 99% confidence 
level (from an F-tcst). The fitted sinusoids also provide 
marginal evidence for a trend towards larger amplitude 
at longer wavelengths (see table ||). Even less certain 
is the possible "eclipse" feature present in the B band 
data at phase^^O.O, which, if excluded, docs lend fur- 
ther support to a trend in the smooth modulation am- 
plitude. However, it is clear that further observations 
from a larger telescope with better seeing are needed to 
confirm these tentative results. 



4 ORIGINS OF THE OPTICAL FLUX 

The detection of a plausibly phased two hour modula- 
tion confirms our identification of the correct counter- 
part, V4580 Sagittarius, in both sets of observations. 
However, we have noted that its X-ray to optical lumi- 
nosity ratio is very small, which implies that it is sur- 
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Table 2: Parameters of sinusoid fits to optical modulations in V4580 Sgr/SAX J1808. 4-3658. 



Lightcurve 


Phase of 


Amplitude (mag) 


Sinusoid vs. constant fit 




minimum light 


Constrained: Period 


Period and Phase 


F-test confidence level 


Whitclight 


0.0084 ±0.0081 


0.063 ±0.003 


0.063 ±0.003 


100.00% 


B band 


0.09 ±0.06 


0.09 ±0.03 


0.08 ±0.03 


99.18% 


B - "eclipse" 


0.14 ±0.06 


0.05 ±0.02 


0.03 ±0.02 


85.23% 


points 










V band 


0.02 ±0.04 


0.08 ±0.02 


0.08 ±0.02 


99.99% 


R band 


0.08 ±0.02 


0.15 ±0.02 


0.13 ±0.03 


99.98% 



prisingly optically bright, especially given the expected 
very faint emission from any secondary capable of fitting 
into the Roche lobe of the 2 hr orbit (see below). The 
possibility that the true quiescent amplitude is much 
larger than in outburst, allowing for another linc-of- 
sight star, responsible for some of the flux we have mea- 
sured, can nevertheless be ruled out by the morphology, 
which is not flat-bottomed but sinusoidal (cf . the case of 
PSR 1957+20, Djorgovski & Evans 1988). We shaU now 
seek to apply some constraints as to the source of the 
quiescent optical emission, given the available informa- 
tion, although we are limited by the large uncertainties 
in the absolute magnitude and colours of V4580 Sgr 
(due to the poorly constrained column). 

Firstly, are we simply seeing the companion star 
and nothing else? Chakrabarty & Morgan (1998) argue 
that the most probable companion is an extremely low 
mass star (< 0.1 ^/q). However, if this lies on the main 
sequence the expected My ~ 14 implies a brightness 
that is 5 magnitudes fainter than our observed limit. On 
the other hand, we might postulate a more massive (0.5- 
0.8 M0) companion, such as a Kl to MO main sequence 
dwarf, which would have the correct magnitude. But 
this is also not tenable. The counterpart is much too 
blue for such a spectral type {{B — V)q < 0.3 cf. a Kl 
star value of 0.86). The inferred inclination would be 
extremely low {i ^ 6°) and therefore improbable, since 
the a priori probability of seeing a system with i or 
less goes as (1 — cosz). Further compelling evidence is 
the lack of ellipsoidal variations. If any modulation is to 
be seen from the light of the companion star alone (as 
is the case for a number of other soft X-ray transients 
in quiescence), then its period should be half that of 
the binary. Hence, we may conclude that no "ordinary" 
lower MS star provides all the optical emission. 

The other possible source of optical emission is a 
combination of: (i) a remnant accretion disc, (ii) a hot- 
ter irradiated accretion disc, (iii) the irradiated face of 
the companion star. In order to make an approximate 
quantitative comparison with our data we have under- 
taken some simplified modelling to obtain predicted 
broad band spectra. For the details of the modelling 
see Chakrabarty (1998), §4.3 (accretion disc), and the 



appendix (heated face of companion star). The resulting 
spectra are presented in figure |^. 

Assuming that an accretion disc is still present, how 
bright will it be solely due to internal viscous energy 
release? Here we make use of the results of the a disc 
prescription, using a nominal but reasonable value of 
q;=0.1, and a canonical 1.4 Mq mass for the neutron 
star. The inner disc radius is set at the co-rotation ra- 
dius assuming a Keplerian disc and the outer at 9/10 of 
the primary's Roche lobe. We may estimate the mass 
transfer rate assuming that the low mass secondary 
continues to fill its Roche lobe (requiring either/both 
an evolved or/and an X-ray irradiation-bloated state), 
and that the process is solely driven by angular mo- 
mentum loss due to gravitational radiation (GR; Kraft, 
Matthews & Greenstein 1962): 

Mgr =i 0.95 X 10"" Afoyr^i (A/2/0.05 Mo)^. 

This is also consistent with the time-averaged value 
based on the energy released during outbursts (e.g. 1998 
April, total fluencc~ 4.3 x 10~^ ergcm"^; Gilfanov et al. 
1998) along with their ^20 month recurrence cycle. The 
spectrum is plotted as the dotted line in the figure. With 
a measured value for ai sin i from pulse timing, the value 
of M2 we choose immediately constrains the inclination 
angle i. For a 0.05 M0 secondary, we find a moderately 
high inclination i = 60° (50% a priori probability of 
detection), and the model curve is consistent with the 
lower limits on the optical fluxes, and has the correct 
blue slope (although M could well be slightly higher 
than we used) . For any higher inclination both the GR- 
driven mass transfer rate decreases (as the square) and 
the projected disc size (as cosi), making such a system 
too faint compared to our observations. On the other 
hand, for more massive secondaries the inclination must 
be lower and hence so will the a priori probability of 
seeing such an inclined system, ignoring any selection 
effects. 

Depending on the system inclination a variety of 
mechanisms can lead to a periodic modulation of the 
optical flux on the binary period. In the case of SAX 
J1808. 4-3658, the lack of an X-ray eclipse requires i < 
82°, whilst we have seen it is more likely to be lower still, 
i < 60°. In this range low amplitude (at most a few per- 
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Figure 5: Optical {BVR) photometry of V4580 Sgr 
(SAX J1808.4-3658) (points with crrorbars), with 
model spectra from various (irradiated) disc/star emis- 
sion components (as labelled). Results for a secondary 
mass of O.OSMQare shown, which immediately deter- 
mines the inclination to be i = 60.3°. The dashed er- 
rorbars indicate the maximum and minimum fluxes ex- 
hibited during the 2 hr modulation cycle, to be com- 
pared with the variation between total from irradiated 
disc (solid line) and total from irradiated disc -I- star 
(long dashed line). No attempt has been made to fit 
the model to the data, instead reasonable values of 
M ~ 10"" Mq yr~i and Lirr ^ 
chosen for illustrative purposes. 



10"^^ ergs ^ have been 



Figure 4: Phase folded hghtcurves [BVR) of V4580 Sgr 
(SAX J18G8. 4-3658), with stars of comparable bright- 
ness plotted beneath each. Fitted sinusoids with con- 
strained period (=1.0) are overplotted on the data for 
V4580 Sgr. For clarity two orbital cycles are shown and 
each curve has been offset by 1 mag. 

cent) sinusoidal modulations are seen in many LMXBs, 
which are presumed to arise due to the varying contri- 
bution from the X-ray heated face of the secondary as 
its aspect changes. Within this framework the measured 
amplitude of modulation enables us to estimate the level 



of irradiating flux required to be Li 



than that of the two quiescent X-ray observations of 
SAX J1808. 4-3658 but probably still consistent, given 
the degree of variability seen in other quiescent sources 
(e.g. Cen X-4; Campana et al. 1997). We note that this 
too argues against lower inclinations, as the amplitude 
will decrease along with the relative contribution of the 
heated face, and a greater degree of irradiation is not 
possible. Interestingly, the trend in amplitude fits that 
implied by the data, decreasing towards the blue, which 
is a consequence of the cooler variable companion star 
component. 



Admittedly, this luminosity is a factor of ten higher 



© 0000 RAS, MNRAS 000, 000-000 



The optical counterpart to SAX J1808. 4-3658: observations in quiescence 7 



5 CONCLUSIONS 

We have undertaken the first time-series CCD photom- 
etry of V4580 Sagittarius, the optical counterpart to 
SAX J1808.4-3658, in quiescence. We find that it has 
dropped to V 21, five magnitudes from the peak of 
the 1998 outburst, although its colours are the same 
within the uncertainties. Furthermore, as during out- 
burst it exhibits a clear 2hr sinusoidal modulation with 
a ^^6% semi-amplitude, and is at a minimum when the 
pulsar lies behind the companion. Our BVR photom- 
etry reveals a possible trend in amplitude, increasing 
with wavelength. Lastly, a comparison of its X-ray to 
optical flux ratios in outburst and quiescence shows 
that if reprocessing of the observed X-ray flux is the 
dominant mechanism in outburst, this cannot apply in 
quiescence, as it is relatively optically over-bright. 

An analysis of these results has led us to conclude 
that: 

• Unlike other soft X-ray transients in quiescence, 
the intrinsic emission of the companion star is not the 
dominant contribution to its optical flux (as evidenced 
by the lack of ellipsoidal variations). 

• If we assume that a low level of mass transfer 
(driven by gravitational radiation losses) continues, 
then the disc emission due to internal viscous energy 
release is the main source of the observed optical flux. 
However, the observed X-ray luminosity requires that 
this mass is not being accreted onto the neutron star. 

• The observed sinusoidal binary modulation is most 
consistent with a geometrical (i.e., aspect change) ori- 
gin, wherein the face of the secondary is irradiated by 
the quiescent X-ray flux, and makes a variable contri- 
bution to the optical flux. 

Clearly, with the existing data we have only been 
able to infer many of the physical details of the sys- 
tem and further observational effort on the quiescent 
counterpart is still needed. 
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APPENDIX A: DERIVING THE 
REPROCESSED OPTICAL FLUX FROM 
THE SECONDARY'S HEATED FACE 

In deriving the optical flux due to the irradiated face of 
the secondary star, a spherical geometry was assumed, 
such that the surface need only be divided into annuli, 
having an angle (^) between their normal and the radia- 
tion incident from the primary. In this way, the method 
used for the accretion disc can be modified for the star. 
Hence, the effective temperature is given by: 

cos^ 

where 77^, is the irradiating flux albedo (taken as 0.3; Mil- 
grom & Salpeter 1976) and d is the separation of the 
annulus from the source. The spectrum is then given 
by assuming each annulus emits as a blackbody at this 
temperature, and summing the contribution from each 
at a given frequency (the shadowing effect of the disc 
has been taken into account). However, for simplicity 
the line of sight is taken as lying in the orbital plane, en- 
abling straightforward calculation of the projected sur- 
face area of each contributing annulus. 
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